GLOBAL EXISTENCE OF SOLUTIONS OF THE LIQUID 
CRYSTAL FLOW FOR THE OSEEN-FRANK MODEL 



Min-Chun Hong and Zhouping Xin 

Abstract. In the first part of this paper, we estabUsh global existence of solutions 
of the liquid crystal (gradient) flow for the well-known Oseen-Frank model. The 
liquid crystal flow is a prototype of equations from the Ericksen-Leslie system in 
the hydrodynamic theory and generalizes the heat flow for harmonic maps into the 
2-sphere. The Ericksen-Leslie system is a system of the Navier-Stokes equations 
coupled with the liquid crystal flow. In the second part of this paper, we also prove 
global existence of solutions of the Ericksen-Leslie system for a general Oseen-Frank 
model in R^. 



1. Introduction 

A liquid crystal is a state of matter intermediate between a crystalline solid and 
a normal isotropic liquid. Research into liquid crystals is an area of a very suc- 
cessful synergy between mathematics and physics. There are a lot of analytical 
and computational issues, which arise in the attempt to study static equilibrium 
configurations. Numerical and experimental analysis has shown that equilibrium 
configurations are expected to have point and line singularities ([K]). Mathemati- 
cally, Hardt, Kinderlehrer and Lin in their fundamental papers [HKLl] and [HKL2] 
proved the existence of an energy minimizer u of the liquid crystal functional and 
showed that a minimizer u is smooth away from a closed set E of $7. Moreover, E 
has Hausdorff dimension strictly less than one. In [AL], Almgren and Lieb did some 
related analysis indicating that the phenomenon is of wider interest. In physical 
theory, an equilibrium configuration corresponds to a critical point, not necessarily 
an energy minimizer, of the liquid crystal energy. Critical points are much harder to 
understand mathematically than minima. From the above result of Hardt, Kinder- 
lehrer and Lin, minimizers cannot have line singularities. Following the work of 
Bethuel-Brezis-Coron on harmonic maps in [BBC], Giaquinta, Modica and Soucek 
[GMS2] found a relaxed energy for the liquid crystal systems, whose minimizers are 
also equilibrium configurations. On the other hand, Giaquinta, Modica and Soucek 
[GMSl] also proved that minimizers of the relaxed energy for harmonic maps are 
smooth away from a 1-dimensional singular set. Further developments on the reg- 
ularity results on harmonic maps were surveyed in [GMS3] . There is an interesting 
open problem to prove that minimizers of the relaxed liquid crystal energy have 
line singularities. The first author in [Ho3] proved partial regularity of minimizers 
of the modified relaxed energy of the liquid crystal energy. However, the partial 
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regularity of minimizcrs of the relaxed energy for liquid crystals is still mysteri- 
ous. In some related studies of liquid crystals, Bauman, Calderer, Liu and Phillis 
[BCCP] studied the Landau-de Gennes free energy used to describe the transition 
between chiral nematic and smectic liquid crystal phase, Lin and Pan [LP] used 
the Landau-de Gennes models to investigate the magnetic field induced instabili- 
ties in liquid crystals, and the existence of infinite many liquid crystal equilibrium 
configurations prescribing the same boundary was obtained in [Ho2] . 

A general description of the static theory of liquid crystals is given by Erickscn in 
[Er] . A liquid crystal is composed of rod like molecules which display orientational 
order, unlike a liquid, but lacking the lattice structure of a solid. The kinematic 
variable in the nematic and cholcstcric phase may be taken to the optic axis, which 
is a unit vector field u in a region C M"^ occupied by the materials. The liquid 
crystal energy for a configuration u e H^{fl; 5^) is given by 

(1.1) E{u;fl)= [ W{u,S7u)dx, 

Jn 

where the Oseen- Prank density W{u, Vu), depending on positive material constants 
ki, k2, ks and ^4, is given by 

W{u, Vu) = fci(div + k2{u ■ curl u)^ -|- ks\u x curl + fc4[tr(Vu)^ - (div uf]. 
Without loss of generality, as in [HKLl] or [GMS3], we rewrite the density 

(1.2) W{u,Vu) = a\S7u\'^ + V{u,Vu), a = mm{ki,k2,k3} > 0, 
where 

V{u, Vu) = (fci — a) (div u)^ + (^2 — a){u- curl u)^ + {k^ — a)\u x curl up. 

A static equilibrium configuration corresponds to an extremal (critical point) of 
the energy functional E in H^{n, S"^). The Euler-Lagrange system for the general 
Oseen-Frank functional (1.1) (see details in Appendix) is: 

(1.3) 

Va [Wpi^ {u, Vu) - u'-uWpi^ {u, Vu)] - W„i (it, Vu) + Wui {u, Vu)u^u' 
+ Wpi^ {u, Vu)Vocu''u'' + Vpi^ {u, Vu)u' Vau' = in O 

for i = 1,2,3, where we adopt the standard summation convention. In a special 
case of ki = k2 = fcs, the system (1.3) becomes the harmonic map equations into 
S"^. However, the equilibrium system associated to the energy functional (1.1) is 
not elliptic for every choice of the constants ki, k2 and k^. 

In the first part of this paper, we investigate the liquid crystal flow for a model 
with the Oseen-Frank density (1.2). For a domain O in E'^ or in M^, a map u{x, t) : 
ft X [0, 00) is a, solution of the liquid crystal flow if u satisfies 

(1.4) 

^ =Va [Wpi^ {u, Vu) - u^u%,^ {u, Vu)] - W„* {u, Vu) 

+ W^l {U, Vu)v}u^ + Wpl^ [U, Vu)^aU^U^ + Vpl^ {U, Vu)vl-VaU^ 
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in Q X [0, oo) for i = 1,2, 3. 

The flow equation (1.4) is a prototype of equations from tlie Erickscn-Leslie 
system in the hydrodynamic theory (of. [Er]). The hquid crystal flow (1.4) also 
generalizes the heat flow for harmonic maps into the 2-sphere. Since the seminal 
work of Eclls-Sampson [ES], many work on the heat flow for harmonic maps have 
been made. In 2 dimensional case, Struwe [Stl] established global existence of the 
weak solution of the harmonic maps flow with initial data, where the solution is 
smooth except for a finite number of singiilaritics. In higher dimensional cases, 
Chen-Struwe [CS] proved global existence of partially regular solutions to the har- 
monic map flow. Since (1.4) is not parabolic, the system of the liquid crystal flow 
is complicated, so the question on global existence for the liquid crystal flow (1.4) 
for the Oseen-Frank model remains unresolved. In this paper, we prove global 
existence of solutions of the liquid crystal flow in 2D. 

We set 

H^{M.'^;S^) := {u: u - b G H^R^iM.^), |w| = 1 a.e. in M^} 

for a constant vector b G S"^. 

Then, one of our main results in this paper is the following global existence for 
this flow in 2D (i.e. w is a constant along a direction in R'^): 

Theorem A. Let uq G if^(IR^;S'^) be a given map. Then there exists a global 
weak solution u{x,t) : x [0, +oo) — S"^ of (1.4) with initial value w(0) = uq 
such that u is smooth in x [0, +oc) except for a finite number of singularities 
{(.t', T|)}/£j^ G X [0, +oo) with an integer K > depending on Uq. Moreover, 
there are two constants £o > and Ro > such that each singular point x\ at the 
time Ti is characterized by the condition 

limsupi? (m (i) , i?K (a;-)) > £q 

for any R> with R < Rq. 

This result can be regarded as an extension of the well-known result of Struwe 
in [Stl] on the heat flow for harmonic maps in dimension two. Since the liquid 
crystal flow is not a parabolic system, the flow (1.4) is more complicated than the 
harmonic map flow. In particular, we can not apply the well-known theory of partial 
difFcrcintial equations directly to prove the local existence for the liquid crystal flow. 
Instead, we consider a family of Ginzburg-Landau approximation flows to prove 
local existence of solutions to (1.4). To prove Theorem A, we need to get a L^- 
estimate of V^tt similarly to one in [Stl]. However, the flow (1.4) is not a parabolic 
system, so we overcome the difiiculties due to the term VaWu^Vpi^ {u, Vw)] by using 
the fact that \u\ = 1 as observed in [Hoi]. 

In the second part of this paper, we investigate the Ericksen-Leslie system with 
the Osccn-Frank density W{u,\7u) in (1.2). In the 1960's, Ericksen [Er] and Leslie 
[Le] established the hydrodynamic theory of liquid crystals independently. The 
Ericksen-Leslie theory describes the dynamic flow of liquid crystals, including the 
velocity vector v and direction vector u of the fluid. Let v = {v^,v'^,v^) be the 
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velocity vector of the fluid and u = {u^ ,u^) the unit direction vector. The 
Ericksen-Leslie system in x [0, oo) is given by (e.g. [LI] and [LLl]) 

(1.5) v\ + {v -vAv' + V^,P= -AV^, {V^.vl'W^H [u, Vu)), 

(1.6) V-u = 0, 

(1.7) 

u\ + {v- V)m* =Va [Wp^ (u, Vu) - u''u%u (u, Vu)] - (u, Vu) 

+ (m, Vu)u^u'^ + VFpI^ (U, Vu) Vaw'w* + Vpfe (u, Vu)u^WcxU^ 

for z = 1, 2, 3, prescribing the boundary condition 

(1.8) v{x,t)=Q, u{x,t) =uo{x), y{x,t) Gdnx{0,oo) 
and with initial data 

(1.9) v{x,0) = vo{x), u{x,0) = uo{x), div t;o = Vx e f2. 

Here v, A are given positive constants, and P is the pressure. 

The system (1.5)-(1.7) is a system of the Navier-Stokes equations coupled with 
the liquid crystal flow (1.4). The study of the Navier-Stokes equations is of great 
interest. Tremendous results on the existence and partial regularity for the Navier- 
Stokes equations have been established (e.g. [Sc], [CKN], [L2], [TX]). In this paper, 
we are only concentrating on the existence of solutions of the Ericksen-Leslie system. 
Since the functional E{u; O) in (1.1) with the constraint l^l = 1 is complicated, one 
considers Ginzburg-Landau functionals 



n 



W{u,Vu) + ^{l-\u\'f 



dx 



for any function u € H^{Cl;M.^). Then, the approximating Ericksen-Leslie system 
is given by 

(1.10) vi + {v ■ Wy -uAv' + V^,P = -AV^,(V^,w*=VFpj(«, V«)), 

(1.11) V-v = 0, 

(1.12) ui + iv Vy = V„ [Wp.^ {u, Vu)] - W^i [u, Vu) + ^u\l - |«|2) 

for i = 1,2,3, prescribing the boundary condition (1.8) and initial condition (1.9). 

In the case of ki = k2 = k^, Lin and Liu [LLl] proved global existence of the the 
classical solution of (1.10)-(1.12) with (1.8)-(1.9) in dimension two and the weak 
solution of the same system in dimension three. Lin and Liu in [LL2] also analyzed 
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the limit of solutions {v^,u^) of (1.10)-(1.12) as e 0, but it is not clear that 
the limiting solution satisfies the original Ericksen-Leslie system (1.5)- (1.7) with 
|u| = 1. Therefore, there is an interesting question to establish the global existence 
of solutions of (1.5)-(1.7) with (1.8)-(1.9). The question for the case of ki — k2 = 
has been answered by the first author in [Ho3] in and Lin-Lin- Wang [LLW] in 
a general case for a domain of independently. The system (1.5)-(1.7) or (1.10)- 
(1.12) for the general Oseen-Frank model is more complicated than the system for 
the case of kx = k2 = since there is no maximum principle for the parabolic 
system (1.12) in the case fci 7^ k2 (see [A]) and the term Wui{u,\7u) in (1.12) will 
cause a trouble to prove global existence for the system. 

In this paper, we will prove global existence of weak solutions to the Ericksen- 
Leslie system (1.5)-(1.7) for a general Oseen-Frank model in M^. More precisely, 
we have 

Theorem B. Let (uo,fo) G H^{R^;S^) x L^(]R^,M^) be given initial data with 
dvi vq = 0. Then, there exists a global weak solution {u,v) : x [0, +00) 

X of (1.5)-(1.7) with initial values (1-9), where the solution {u,v) is smooth 
in X ((0, +oo)\{T;}^-|^) for a finite number of times • Moreover, there are 

two constants £0 > and Ro > such that each singular point (a;-,T;) e S x {T;} 
is characterized by the condition 



for any R> with R < Rq. 

The main idea to prove Theorem B is to combine the idea in [Ho3] with the 
proofs of Theorem A. The first key step is to prove local existence of solutions of 
the system (1.5)-(1.7) by considering the approximation system (1.10)-(1.12). To 
prove global existence of solutions to (1.5)- (1.7), one of key steps is to get a L^- 
estimate of V^u and Wv in M? x [0,r] under a small energy condition as in [Stl]. 
To show the regularity of the weak solution {u,v) of (1.5)-(1.7) in x (0,T), we 
establish a local energy inequality under the small energy condition, which was 
first used by Struwe in [St2] for the iJ-system flow. Finally, we prove regularity 
of solutions by controlling L^-estimate of V^u and Vv in M.'^ for t E (0,T). Since 
(1.7) is not a parabolic system, the proof of Theorem B is more difficult than one 
for the case of fci = ^2 = k^ in [Ho3]. We overcome a number of difficulties on 
the regularity and uniqueness for the systems by employing the invariance of the 
density (1.2) after a rotation. 

The rest of the paper is organized as follows. In Section 2, we prove the global 
existence for the liquid crystal flow in 2D. Some global estimates for (1.5)-(1.7) arc 
established in Section 3. Then, we complete a proof of Theorem B in Section 4. 
Finally, the regularity issue for the systems is dealt in Section 5. 

2. Existence of pcirtial regulcir solutions of the liquid crystal flow 

In this section, we consider the flow (1.4) in M?. For simplicity of notations, u is 
assumed to be a constant along ars-direction in M^; i.e. -S^ = 0. 
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For any two positive constants r and T with t <T,we define 
V{t, T) := {u •.M.^x [t, T] — >• 5^ I u is measureable and satisfies 

ess sup / \S7u{-,t)fdx+ [ [ \S7'^uf + \dtuf dxdt < 00} . 
T<t<T Jr2 Jt Jr^ 

Lemma 1. Let u € V{Q,T) he a solution of the system (1-4) with initial value 
uo e 52). Then, for any ti e [0,T] 

(2.1) [ \dtu\'^dxdt + E{u{ti)) <E{uo). 

JK2x(o,ti) 

Moreover, for all t G [0,T], Xq e and R> 0, it holds that 

(2.2) / W{u{x,t),Vu{x,t))dx 

■>Br{xo) 

< I W{uo{x),Vuo{x))dx + C^ / \Vuo\^dx, 

where C is a constant. 

Proof Multiplying (1.4) by ^ yields 

/ \ — \^dx = - Wp,{u,Vu)—Vau'dx- W^^{u,Vu)—dx. 
Jr2 at at at 

This implies 

Jm? (Jt at J'g2 

(2.1) follows from integrating the above identity. 

Let (j) e C^{B'2r{xo)) be a cut-off function satisfying < < 1, |V0| < C/R 
and (f)=lon Buixo). Multiplying (1.4) by and then using Young's inequality 

yields 

/ \^\^^^dx + ^ [ W{u{x,t),Vu{x,t))cP^dx<C [ |^||V«||<^V<^|da; 

'■dx + C I W{u,\/u)\\I<j)\'^ dx. 

Then, (2.2) foUows from using (2.1) and integrating the above inequality. □ 

It follows from [Stl] that 

Lemma 2. There are constants C and Rq such that for any u S V{0,T) and any 
R e (0, Rq], we have 

/ \Vu\'^dxdt<C esssup / \\7u{-,t)fdx 

■{[ \S/^u\'^dxdt + R-^ [ \S/u\'^dxdt). 

Jv? X [0,T] Jw X [0,T] 



< 1 / 1^12^2. 

- 2 7r J dt 
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Lemma 3. Let u G V{0, T) he a solution of (I.4) with initial smooth value uq G ■ 
Then there are constants e\ and Rq > such that if 



essE 

0<t<T,'. 



isup / \Vu{-,t)\'^ dx < si 

T,xeM2 Jbr{x) 

for any R e (0, Rq], then 

(2.3) / \y^u\^dxdt<CE{uQ){l+TR-'^), 

Jv? X [0,T] 

(2.4) / \Vu\^ dx dt < CeiE{uo) (1 + Ti?-^). 

Jr^ X [0,T) 

Proof. Multiplying (1.4) by Au* yields 

= / Va [Wpi {u, Vu) - u'^u'Vpk {u, Vu)] A u' dx 

- [ W^.{u,Vu){Au' -u'^u' Au'')dx+ [ Wp<= (w, Vu)VauV A uMa; 

+ / Vpk {u, Vu)u''Vau' A da; := /i + /2 + /a + /4- 

Note that the terms I2 and of the above identity can be controlled by C|VmP| Am|. 
It suffices to estimate terms Ii and I4. Since = 1, — A = |VMp. We note 

Va[u''u'Vpk{u,Vu)] =\7au''u'^Vpk{u,\7u) + u''\/au'Vpk{u,Vu) 

+ u''u'\/aVpk (m, Vm). 

Integration by parts twice yields 

h+h= [Wp. (w, Vw)] Wl/^u' dx+ Vau'^u'Vpk {u, Vu) A u' dx 

- I u^Voypk{u,Vu)\Vu\'^ dx. 



Note 



VoiVpk {u, Vw) = Vpk pi {u, Vu)V^„w' + Vpkui {u, Vu) Vcw' 



and 



VpWp^^ {u, Vu) = Wp^p. (u, Vw)V2^w^' + (w, Vu)V0U^. 

This implies 

(2.5) ^/ \Vu\Ux+ [ W^. ,{u,Vu)Vl^uW^^puUx 

<C [ |Vu|2(|Vw|2 + |V2w|)rfx. 



Since W(u,p) is convex in p, it satisfies the ellipticity 

Wp,^pi^iu,Wu)Wlpu'w'^0U^ > alW'uf 

for the constant a > 0. Then, choosing ei > to be sufficiently small and applying 
Lemma 2 lead to (2.3) and (2.4). □ 
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Lemma 4. Let u G V{0,T) be a solution of (1-4) with initial value uq G H^. 
Assume that 



ess sup / \Vu{-,t)f dx < ei 

0<t<T,xeM2 JBr{x) 

for any R € (0, Rq] . Let r e (0, T] be any constant. Then it holds for all t £ [r, T] , 

(2.6) / \V^u{x,t)\^ dx <Co, 

with a uniform constant depending only on t, T, Rq, and E{uo). 

Proof. The proof is similar to [Stl; Lemma 3.10]. Using a proper cut-off function 
if necessary, we assume in the following proof that / IdtV u\'^{-,t),dx is finite. 

Differentiate (1.4) with respect to t, multiply the resulting identity by 9tu', and 
then integrate to obtain 

f \dtu\^dx+ [ W^^.{u,\Iu)Vadtu'VpdtuUx 
I at Jr2 Jr2 ^'^^li 

<C [ {\dt |Vu|^ + |Vu| \dt u\ I V^t u\] dx 
+ [ dtu' Va [u' Vpk {u, Vdt u)] dx. 

Due to the convexity of W{u,p) in p, there exists a positive constant a > such 
that 

/ W.j{u,Vu)Vadtu'Vf3dtU^dx<a \Vdt u{x,t)\'^ dx. 
Since |u| = 1, so dt = 0. And hence, 

[ dtu'Va {u' Vpk {u, Vdt u)) dx<C [ \dt u\ \Vu\ | VS* u\ dx. 

It follows from these and Cauchy's inequality that 

(2.7) ~[ \dtu{x,t)\^dx+^ [ \Vdtu{x,t)\^dx<C I \dtu\^\Vu\Ux. 

I dt Jj^2 2 J^2 Jj^2 



Note that 



C [ \dtu\^ \Vu\^dx 
<c(^J \dtu\Uxy (^j iVul^da;^ 

<c(^j \dtu\^{x,t)d:^' (^j \dtVu\'' {x,t)^' (^j IVwI^rf^y 

<^/ \Vdtu{x,t)\'^dx+{c [ \\Iu{x,t)\'^ dx] f \dtu{x,t)\'^ dx. 

4 Jr2 V iR2 / 
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This, together with (2.7), yields that for all t e (0,T], 

(2.8) 4/ \d,u{;t)\^dx+^ f \mu{;t)\^dx 

<(C f \Vu{-,t)\Ux] [ \dtU{;t)\^dt. 

It follows from (2.8), (2.4), Lemma 1, and Gronwall's inequality that for any < 

s <t<T, 

f |atn(-,t)|2(it< (e^/a/«2|v«(.,;)rd.d/\ j \d^u{;s)\'' dx 

Combining this with (2.1) shows that for any fixed < r < T, there exists a 
constant C such that 



esssup / \dtu{-,t)\'^dt<CT-'^ E{uo)e'^^^' 

r<t<T Jr2 

with a uniform constant C. On the other hand, using (2.5), integration by parts 
yields that for any t G [r, T] , 

[ \V^u{;t)fdx<C [ \Vu{-,t)\^dx + C [ \dtu{-,t)fdx 

Jr2 JR2 Jr2 

<Csi[ \V^u{;t)fdx + ^E{uo) + c[ \dtU{;t)\^dx. 

Jr^ -Kq Jr2 

Combining this with (2.9) shows that for suitably small ei, the desired estimate 
(2.6) holds with 

(2.10) Co = CE{uo) + e^'^ B(«o)(i+Tij-^)^ _ 

By the well-known Gagliardo-Nirenberg-Sobolev inequality, we have for any x € 
\u{x,ti)-u{x,t2)\ < C||u(a;,ii) -u(a;,t2)||^4%^(^))||u(a;,ti) -M(a;,t2)||i''2*B^(^)) 

/ T \ 

<C{ sup \\VM-M%\R2. + l)\ti-h\'/'l [ [ \dtu\^dxdt] 

T<t<T ' \Jo Jr^ J 



< C\h - h 



11/8 



It follows from (2.6) and Sobolev embedding theorem that u{x,t) is Holder con- 
tinuous in X uniformly for t G [t, T]. Then we get that u is Holder continuous in 
C^/^(K2 X [t, T]) for any T <Ti. Due to Proposition 14 in Appendix, u is in 
Hence, u is regular in (0,Ti). □ 
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Remark 5. Let u G V{Q,T) he a solution of (1-4) with initial value uq € H^. 
Assume that there are constants ei and Rq > such that 

esssup / \Vu{-,t)\'^ dx < si 

0<t<T,£cGR2 J Bb.{x) 

for any Rg {0, Rq]. Then, for any t € [0, T] and R< Rq, we have 
\V\{x,t)fdx <Ci = Ci{\\uo\\m,Co). 



J 



Theorem 6. (Local existence) For a map uq G i/^(M-^,5^), there is a solution 
u G V{Q,t\) of (1-4) with initial value uq for some ti > 0. 

Proof. For any map uq e {M."^ , S"^) , it can be approximated by a sequence of 
smooth maps in H^{R'^,S'^). Without loss of generality, we assume that uq G 
H^{M.'^,S'^) is smooth. The liquid crystal flow is not a parabolic system, so one 
can not apply the well-known local existence theory. Instead, we prove the local 
existence by an approximation of the Ginzburg-Landau flow in the following: 

(2.11) -dt = ^" [Wp^Ju„Vu,)] - W^.{u„Vu,) + -<(1 - \u,\^) 

with initial value Ho G _ff^(R^,S'^) and uo G C°°. Applying the standard local 
existence theory of quasi- linear parabolic systems (of. [Ei] or [Am]), there is a local 
regular solution of (2.11) with initial value ^£(0) = uq- 

For simplicity of notations, we deflne 
y(r, T) = |u : X [r, T] ^ IR^|w is measurable and satisfies 

esssup / \\/u{-,t)\^dx+ [ I {{V"^ + \dtu\^)dxdt <oo\, 
T<t<T jR2 Jt Jm2 J 

< T < T < +00, 
e,{u) = W{u, Vu) + 7:^(1 - |w|')', E,{u) = [ e,{u) dx. 

Taking inner product of (2.11) with dt u^, one can obtain that for any s > in the 
maximal interval of existence, 



/ 



(2.12) / \dtUe\^dxdt + E,{Ue{s)) <E{uo). 

<(0,s) 



Moreover, repeating similar arguments in Lemma 3 (see below (2.29) below) yields 
that the solution belongs to V"(0, T^) for a maximum time and hence is regular 
X [0, Tg). The maximum time Tg is characterized in the following: For a singular 
point xo at T^, there are eq and iio > such that 

limsup / \\/ue{-,t)\'^ dx > So > 



THE LIQUID CRYSTAL FLOW 11 

for any positive R< Rq. 

Next, wc will show that there is a uniform lower bound time ti > such that 
Te > ti and is bounded in V{0,ti) uniformly in e. 

A similar argument as in Lemma 1 shows 

/ ee{ue{x,t)) dx < / es{uo{x)) dx + C ^ / \Wuo\^ dx 

JBr{xo) JB2r{xo) ^ Jv? 

for t < Te. 

It follows from this inequality that for suitably small e\ and -Ro, there is a time 
ti uniform in £ with t\ < Tg such that 

(2.13) sup / ee{ue{x,t))dx < £i 

0<t<ti Jbr{xo) 

for R < Rq and thus is smooth for [0, ti] for all e > 0. Next, we claim that for 
0<t<ti 

1 3 

- < \ue{x,t)\ < - for all x e R^. 

To verify this claim, we re-scale the solution by u{x,t) = Ue{ex,e'^t). Then u 
satisfies 

(2.14) = ^« - + - 1^1') 

with initial value Uo{sx). Let r be the maximal time in [0, |i] such that (2.14) 
holds, i.e., 

(2.15) l<\u{x,t)\<^ 

for any {x,t) gM? x [0,t]. Note that in this case, the basic energy inequality (2.12) 
becomes 

(2.16) 

/ \dtufdxdt+ [ {W{u,Vu){s) + l-{l-\u{s)\'^f)dx < E{uo) 



([0, 

for all s e 



and the condition (2.13) turns into 

(2.17) esssup I ( |Vw(-,s)|^ + ^(1 - |u|^)^) rfa; < £i. 



Q<S<^f,X&? J Br {x) 

for R < Ro. 
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Multiplying (2.14) by Au and integrating over lead to 

(2.18) ~[ \\/u\'^dx+ [ Va (u, Vu)] A wMx 

- / Wj,i{u,Vu)Au'dx+ [ u\l-\u\'^) Au'dx = 0. 



Note that 



/ u'(l- lup) AwMx = - / |Vup(l- |w|2)dx + 2 /" \V\u\''\'' dx. 
Jk? Jr^ Jr^ 

Then, combining the above identity with (2.18) yields that for any s,t G [0, r] with 
s<t, 

(2.19) 

I f \Wu{;t)fdx+ f I W^,^{u,^u)^l^u'W'^^uUxdt 



< 



cj \Vu{;s)\^dx + C [ I \yu\'' + ri[{l-\u\^f + \^^u\^]dxdt 

Jr? Js Jr^ 



for a sufficiently small 77 > to be chosen. 

On the other hand, it follows from (2.14) and (2.15) that 

(2.20) / {l-\u\^fdxdt<c[ i\dtu\^ + 1^111"^ + \\7^u\)dxdt. 

jR^x[s,t] JR^x[s,t] 

Combining Lemma 2 with (2.17) shows that 
(2.21) 

/ \Vu\^dxdt<CiSii [ \\/'^u\'^dxdt + ^ [ \yu\^dxdt]. 

jR^x[s,t] \JR^x[s,t] -Rq JR^xls,t] J 

As a consequence of (2.19)-(2.21), (2.16), and suitable choices of r] and ei, one 
can get that u € Vr, and for any < s < t < t, 

(2.22) / [\V^u\ + {l-\uff]dxdt<CE{uo){l + e'^{t-s)Ro'^), 

jR^x[s,t\ 

(2.23) / \'S/u{x,t)\^dxdt <CsiE{uo){l + s^{t- s)Ro'^). 

jR^x[s,t] 

By a similar argument as in the proof of Lemma 4, one can derive from (2.14) 
that there exists a positive uniform constant a such that 

(2.24) 

~f \dtuix,t)\Ux + a [ \dtyuix,t)fdx+l f \dt{\uix,t)\^)\^ dx 

<C I \dtu{x,t)\^\Vu{x,t)\^dx+ I \dtu{x,t)\^ {l-\u{x,t)\^)dx 
Jr^ Jr^ 
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Note that 

C I \dtu{x,t)f \Vu{x,t)fdx 

<c{^j \dtu{x,t)\^dxy (^J \Vu{x,t)\^dxy 

1 1 , - 

<c([ \dtu{x,t)\'^dx){ f \dtWu{x,t)\'^dxy ( f \Vu{x,t)fdx 
<^ f \dtVu{x,t)\^dx + c( f \Vu{x,t)\'^dx) f \dtu{x,t)\'^ dx. 



d_ 

dt , 



Hence, 

/ \dtu{x,t)fdx + a [ \dt'^u{x,t)f dx + [ \dt{\u{x,t)\^)f dx 

Jr2 Jr2 Jr2 

c[ \Vu{x,t)\'^dx) f \dtu{x,t)\^dx + 2 f {dtu{x,t)\^{l - \u{x,t)\'^)d2 

JR2 / JR2 J-R2 

which yields immediately that for any 0<f<rG(0,|y], 

/ \dtu{x,t)\^ dx 
< e^/o/82 \Vu{x,t)\Uxdt I u(x,0)|2 dx 

JR2 

It follows from this, (2.23), Uo G H^, and (2.14) that 
(2.25) / \dtu(,x,t)\^dx<Ct=Cx{E{uo),ei,tu\\uo\\Hi,Ro) 



with a positive constant Ci independent r e (0, |j-) given by 

eii5(«o)(l+|^) 



(2.26) Ci=(7(|Kllff2)e 

Using (2.18), an integration by parts yields implies that for all t G (0,r], 
/ \V^u{-,t)fdx<C [ \Vu{x,t)\^dx+C [ {l-\u\^fdx+C [ \dtu{x,t)\^ dx 

JR2 Jm.2 J]^2 Jr2 

Due to Lemma 2, and (2.17), one has 

c[ \Vu{x,t)\'^ dx < ei [ \V^u{x,t)\^dx + ^^^E{uo). 
Jr2 Jm.2 Rq 
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Thus one can get that for all t G (0,t). 
(2.27) / \y^u{x,t)\''dx<CE{uo)il + ^) + CCi. 

By the Sobolcv embedding theorem, u is /3-H61der continuous in x uniformly in 
all t € [0,t] with /3 < 1. Repeating the similar analysis as in the proof of Lemma 4 
and using Proposition 13 in Appendix, we get u G C^'» on x {0,t). If there is 
a xi G such that cither \u{xi,t)\ < i or |{((xi,f)| > |. By the uniform Holder 
continuity of u, there exists a constant C2 with the property that < and 

{l-\u{x,t)\^f>\, xGB^ix,). 

Hence, 

{l-\u{x,t)\Ydx>^\B^^{0)\>2e^. 

which contradicts to (2.17) for suitably small £1. Here we have used the fact that 
C2 depends only the upper bound of Ci , which may be chosen to be independent 
of £1 by the choice of ti. This implies that 5 < |{t(x,t)| < | for all t G [0,t]. By 
the continuity of m at r and the maximal choice of r, r must be the value |y. This 
shows that (2.14) holds for all t e [0,ti]. 

Next, it follows from (2.12) and (2.22)-(2.23) that Ug are uniformly bounded in 

V{0,ti) for aU e and 



(2.28) 




(2.29) r/ \\^'Mx,t)\' + ^{l-\ue{x,t)ff]dxdt<CE{uo){l + ^), 
Jo JR^ L ^ J -"o 

(2.30) /Y \'^Mx,t)\'^dxdt<CeiE{uo){l + -^). 
Jo Jkp- -^0 

Letting £ — )■ 0, we can prove local existence of a solution of (1.4) in ^(0, ti). □ 

Now we complete the proof of Theorem A. 

Proof of Theorem A. By Theorem 7, there is a local solution u on [0,ii) for some 
ti > 0. By Lemma 3 and Lemma 4, the solution can be extended to [0,Ti) for a 
maximal time Ti > such that there is a singular set S at Ti. Each singularity 
€ S at Ti is characterized by the condition 

limsupii^ (u (t) , Bn {x])) > £0 

t/'Ti 

for any R > with R < Rn- It is easy to sec the solution u £ F is regular for all 
t e (0,ri). By Lemma 1, we can show that the singular set S and the singular 
times are finite (Sec [Stl]). Theorem A is thus proved. □ 

Remark. There is an open problem to prove the uniqueness of the weak solutions. 
But, we can prove the uniqueness of smooth solutions (see Lemma 11 below in 
Section 3). 
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3. Global existence for the Ericksen-Leslie system 

In this section, we derive a-priori estimates for solutions to the Ericksen-Leslie 
system (1.5)-(1.7). Without loss of generality, we assume that f = X = 1 in (1.5). 

For the case Q. = R^ we still consider (1.5)-(1.7) in by taking ^ = 0, ^ = 

0. In this case, V-f; = |f^ + |f^=Oin (1.6) is well-defined. 

For two positive constants r and T with r < T, we denote 

V{t, T) := {u : X [r, T] 5^ | u is measureable and satisfies 

esssup / \Vu{-,t)\'^dx+ [ [ u\'^ + \dtuf dx dt < oo} 

T<t<T Jr2 Jt Jr2 

and 

H{t, T) := {v iR'^x [r, T] ->■ | v is measureable and satisfies 

esssup / \v{-,t)\'^ dx + / / \\7 V \^ dxdt < oo}. 

r<t<T Jr2 Jj. Jr2 

For each pair {u, v), define 

e{u,v) = W{u,\7u) + —\v\'^ , E{u,v) = [ e{u,v)dx. 



M2 



Lemma 7. Let {u,v) e V{0,T) x H{0,T) be a solution of (1.5)-(1.7) with initial 
values uo e H^{M.'^; S'^) and vq G ^^(M^;]^^). Then for t e (0,T], 

(3.1) / e{u{-,t),v{-,t))dx+ I [ {\ut + {v ■V)u\^ + \Vv\'^)dxdt 

Jm? Jo Jr^ 

= / e{uo,vo)dx. 

Proof. Multiplying (1.5) by v and using (1.6), one gets 

(3-2) HI \v\'^dx+ [ \S/v\'^dx= [ Vjv''^iu''Wpk{u,Vu)dx. 

Multiplying (1.7) by Ut + {v ■ V)u yields 
(3.3) 

/ {u\ + {v ■ V)w') (Va [Wp. {u, Vu) - u''u%k {u, Vw)] ) dx 
+ / {ul + {vV)u'){-W^i{u,S7u) + W^k{u,Vu)u''u'')dx 

iR2 

+ / {ui + {v- V)m') ( + Wpk {u, Vu) Vau'^u* + Vpk {u, Vu)u''VaU*) dx 

JR2 " " 

= / \ut + {v- V)u|^ dx. 
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Note that |u|^ = 1 implies 

u^dtu^ = 0, u^^^y = 0. 

Integration by parts yields 
(3.4) 

/ ui (Va [Wpi {u, Vu) - u''u''Vpk {u, Vu)] - Wui [u, Vu) + Vpk {u, Wu)u'"Vau'') dx 
= - Vaul [Wpi {u, Vu)] - ulWu.{u, Vu) dx = -— / W{u, Vu) dx. 

Using (1.6) and integrating by parts, we get 
(3.5) 

/ {V ■ V)U* ( Va [Wpi {u, Vu) - u''u%k {u, Vu)] - W^. {u, Vu) 
+ Vpk {u, Vu)u'^VaU^) dx 
= - / Vav'^VkU^Wpi {u, Vu) + [Vfe Vaw'Wpi {u, Vu) + VfeU W„i {u, Vu)] dx 

= - VaV^Vku'-Wp^ (u, Vu) dx. 
It follows from (3.3)-(3.5) that 



(3.6) 

d_ 
'dt 



I W{u,Vu)dx+ \ut + {v ■ V)u\'^ dx = - Vav''Vku''Wpi {u,Vu) dx. 
Jr^ Jr^ Jr^ " 



Therefore, (3.1) follows from integrating (3.2) and (3.6) in t. □ 

By the same proof as in [Stl; Lemma 3.1], there exists a constant Ci such that 
for any / e H{0, T) and any i? > 0, it holds that 

(3.7) 

/ \f\*dxdt<Ci CSS sup [ \f{-,t)\^dx 

7r2 X [0,T] 0<t<T,xeR^ JBaix) 

■ I [ IVffdxdt + R-"^ [ \ffdxdt] . 

\Jm^xIO,T] Jr^x[0,T] J 

Then, we have 

Lemma 8. Let {u,v) & V{0,T) x H{0,T) be a solution of (1.5)-(1.7) with initial 
values uq G and vo € L^. Then there are constants e\ and i?o > such that if 

ess sup / e{u{-,t),v{-,t)) dx < El 

0<t<T,xeR'^ JBr{x) 
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for any R G (0, Rq], then 

(3.8) / IV^ul^ + l'^vl^dxdt^Cil+TR-^) [ e{uo,vo)dx, 

Jr^ X [0,T1 Jm^ 



(3.9) 



/ {\Vuf + \v\'^)dxdt<Cei{l + TR-^)[ e{uo,vo)dx 

Jm.^ X [0,T] Jr^ 



Proof. Multiplying Au* with (1.7) yields 
/R2 dt 



/ Va [Wp. {u, Vu) - u^u%u (u, Vu)] A w' 

- / W„i (u, Vu) A uMx + / W„fc (u, Vu) A dx 

+ Wpk{u,yu)Wau''u^ Au* dx + Vpk{u,Vu)u'^VaU^ Au'' dx. 



As in the proof of Lemma 3, one can derive 

d 



■ / \Wu\' dx+ f W^, . {u, Vu)Vlf,u'v'^^u^ dx 



dt 

<C f {\Vu\^ + \v\^){\Vu\^ + \V^u\)dx 
b 

< - 

- 4 



/ \V^ufdx + cf {\Vu\^ + \vf)dx. 

JR2 JR2 



/R2 

Applying (3.7) and Lemma 2 again shows 



/ \Vu\^ + \v\'^dxdt 

JR2 X [0,T] 

<Ci£i / |V^up + |V?;pdxdt + Cieii?-2 /" |Vup + dccdt. 

^M^ X [0,T] Jk.^ X [0,T] 

Then (3.8) and (3.9) follow by choosing £1 = 4^- □ 

Lemma 9. Let {u,v) be a solution of (1.5)-(1.7) with initial values {uq,vq) with 
u G F(0, T) and v e H{0, T). Assume that there exist constants ei > and i?o > 
such that 

sup / \Wu{x,t)\'^ + \v{-,t)\^ dx < ex. 

Then for all t € [0,T], xq € IR and R < Rq, it holds that 
(3.10) 

/ e{u{-,t),v{-,t))dx+ [ [ {\Vv\^ + ]-\dtU + v-Vu\^)dxdt 

Jbr{xo) JoJbr{xo) 2 

</ e{uQ,VQ)dx + C2^—{1 + e{uo,vo)dx, 
Jb2r{xo) Jm^ 
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where C2 is a uniform positive constant. 



Proof. Let (f) e Cq°{B2r{xo)) be a cut-off function with = 1 on Bii{xo) and 
|V.^| < IW^cpl < ^ for all R < Rq. 

Multiplying (1.5) by vcly^ and integrating show 

/ vt ■ vcj? + {v ■ V)v ■ vcf)^ - Av ■ v(f? + VP • vcf)^ dx 

= I W^^u^Wpu{u,Vu)V^^v'(t>''dx+ I V^^u^Wpu{u,VuyV^^(f^dx. 
Integrating by parts yields 

[ {vfV(t>^ - Av ■v(t>^)dx=]-^ [ \v\'^(t)^dx+ [ Vv ■V{v(j)^)dx 

= \t I \vf<f>'^dx+ [ \Vv\^(j)'^dx- ( \v\^{\V(l)\^ + (j)A(f))dx. 

2 dt Jr2 7R2 J^2 

Integrating by parts and using (1.6) give 

/ V^^Pv'cf^ dx = -2 I Pv'(t)V^,(i>dx 



and 



\l v'V,,{\v\')<l>' = - [ v'^lvl'd^V^.cl^dx. 



Hence, 
(3.11) 



ld_ 
2di 



f \vm''dx+ f \Wv\''(t>'^ dx 
Jv? J 

{\v\^ +2P+\Vu\^)v'(t>V^,(l)dx+ [ \v\^{\V(l)\'^ +(j)A(j))dx 

V^,«'=Wp.(w,Vu)V^^vV2da;+ j \/^^u''W^^{u,VuyV^.<f'dx. 
Ik? ' Jk? ' 

Multiplying (1.7) by {u\ + {v ■ V)u*)(^^ and using \u\ = 1 lead to 
/ \ut + {v ■ V)up(/>^ dx 

= / {ui + v'\/iu')\/a[Wp. {u, Vu) - u'^u'Vp^ {u, Vu)](f^ dx 

+ / {ui + v^Viu'){-W^i{u,Vu) + Vpk{u,Vu)u''\/au')<l>'^dx. 
Integration by parts yields 

/ wjVa [Wpi {u, Vu) - u^u^Vpk {u, Vu)](j)'^ dx 

+ / uj(- W„i {u, Vu) + Vpk {u, Vu)u''Vau')^^ dx 

= - / [VaulWpi (u, Vw) + WjW^i (w, Vu)]<p'^ dx- ulWpi {u, Vu) Va(?i^ dx 

JR2 " JR2 " 

= [ W{u,Vu)fdx-2 [ vy^p^{u,^u)<|)V„(|)dx. 

dt Jk2 J^2 " 
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Integrating by parts twice and using (1.6), we obtain 

/ v^Viu'ValWpi {u, Vw) - u''u%k {u, Vu)](f)^ dx 

= - I {VaV^Viu' + v^ViVau')Wpi {u, Vu)(f>'^ dx 

- [ {v^S7iu')Wpi {u, Wu)Voc4>'^ dx- j v^Vlu'W^^ {u, SIu)(t>^ dx 

= - I Vc-y' ViuWpi (u, Vu)(j)'^ dx- [ {v^Viu')Wpi {u, Vu)Va^'^ dx 

i/W{u,S7u)(l)yi(j)dx. 



Combing above three identities yields 
(3.12) 

^ [ W{u, Vu)(/>2 dx+ [ \ut + {v ■ V)u\^<l)^ dx 

= - / \7^v''Vku'Wpi^ {u, Wu)(j)'^ dx- f {u\ + v''Wku')Wpi {u, Vu)^^^'^ dx 

+ 2 1 v''W{u,Vu)4>Vk(t>dx 
< - / V^u'^VfeuWpi^ {u, Vu)(j)'^ dx+]- ( \ut + {v- V)up?!>2 

+ C [ \\/u\^\\/<j)\^ dx + 2 [ v''W{u,\/u)<j)\/k(l)dx. 

Integrating (3.11) and (3.12) in t on [0, s] leads to 
(3.13) 

e{u{-,s),v{-,s))(p^dx+ / / {\S/v\^ + l-\ut + {v ■ V)u\^)(l)^ dxdt 

Jo JR2 2 



< 



+ 2 
+ C 



[ e{uo,vo)<l>'^dx+ [ [ {\vf + \Wuf + 2Py(l)V^,(j)dxdt 
Jr^ Jo 

/ / {v ■V)u''Wpk{u,Vu)<pVa:,'pdxdt + 2 [ [ v^W{u,Vu)(l)Vi<pdx 

Jo JR2 * Jo JR2 

r / {\v\' + \Vu\'){m' + \cf>\\A^\)dxdt. 

Jo 

This, together with (3.1), shows immediately that 
(3.14) 

/ {\v{; S)|2 + \Vu{; s)\^) dx + T [ {\Vv\^ + hu^ + {v V)^^^^ dxdt 

Jbr(xo) Jo Jv? ^ 

<! (|uo|' + |Vwo|^)rfa; + C / j {\v\'' + \Vu\'' + \P\)\v\\(l)\\V(l)\dxdt 
Jb2r{xo) Jo JR2 

' f {\vo\^ + \Vuo\^)dx. 

JR2 
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It follows from Holder inequality, (3.1) and (3.9) that 

(3.15) / / {\v\^ + \Vu\^)\v\\(f>\\V(l>\dxdt 

Jo JR^ 

<C r f ilvf + \VufMdxdt 

Jo JR2 -K 

< C (^j^j (|t;|4 + \Vu\'')dxdt^ ' (^^y ^^dxdt 
<C£^—[l + ^j J ^e{uo,vo)dx. 

Similarly, 

(3.16) ^ 

/ / \P\\v\\(l)\\W(l)\dxdt <C f f \p\^-^dxdt 
Jo Js.^ Jo Jm^ R 



<C'^ (^j ^e{uo,vo)dx^ ' ^ J ^ \P\^dxd?j 



for R < Rq. 

On the other hand, it follows from the relation that 

AP = -Vx,x, [Va.,M'=Wp^. {u, Vm) + v^v^ on M2 X (0, T], 
due to (1.5), and the Calderon-Zygmund estimate (cf. [CKN]) that 

/ / \Pfdxdt<c[ [ {\Wuf + \vf)dxdt<Cei(l + ^)f e{uo,vo)dx. 
Jo Jr^ Jo Jr^ ^ R ■' Jr^ 

This, together with (3.16), yields 

(3.17) / / \P\\v\\<p\\W(P\dxdt<C^(l + ^Y ej [ e(uo,vo)dx. 

Jo Jr^ ri \ tt / Jr2 

The desired estimate (3.10) now follows from (3.14), (3.15) and (3.17). □ 

Lemma 10. Let u G F(0,T) and v G H{0,T) be a solution of (1.5)-(1.7) with 
initial value {uo,vo) G Hl{M?, S^) x L'^{M?,M.^) and div vo =0. Assume that there 
are constants e\ and Ro > such that 



sup / \Vu{-,t)f + \v{-,t)\'^ dx < ei 

r,2:GR2 Jbr(x) 



ess 

0<t<T,2:GR2 JBr{x) 

for any R G (0, i?o]. Let r be any positive constant. Then, for t G [r, T], it holds 

that 

(3.18) / \V^u{x,t)\^ + \Vv{x,t)\^dx<CT-^{l + TR-^). 

Jr^ 
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Moreover, u and v are regular for all t £ (0,T). 

Proof. Note that, in a-priority, J^^ I A?;p and J^2 |V^up might not be finite. How- 
ever, by a standard cut-off argument, we can assume that | At;|^ and jV^up 
are finite without loss of generality in the following proof. 

Multiplying (1.5) by Ai>* and integrating by parts, we obtain 
(3.19) 

i-^ / \Vv\'dx+ [ \Av\'dx 

A at Jr2 

= / ivVv')Av'dx+ [ '^j[Viu''Wpk{u,\7u)]Av'dx 
Jb2 7r2 ^ 

] I \Av\^dx + C [ \vVv\^dx + C [ (|V2u|2 + |Vw|^)|Vwpda;. 

4 Jr2 iK2 Jr2 



< 



Differentiating (1.7) in xp, multiplying the above equation by A and then 
integrating by parts, one can obtain 

(3.20) 

-~ [ \Au\^dx+ [ [{Vi3V ■ + {v ■ V) V/3w'] A dx 
2 dt Jr2 J-gp. 

= I [V/sVa \Wpi (u, Vu) - vJ'u'Vpk (m, Vu)] - V^W^i (w, Vu)] A u' dx 

+ I Vfi[Wuk{u,Vu)u''u' + Wp>.{u,Vu)Vau''u']Vi3 Au' dx 

+ I [^p" (■"' Vu)m'' VaW'] V/3 A dx. 

The first term on the righthand side of (3.20) is a bit more complicated. Since 
W{u,p) is quadratic in p, we have 

(3.21) 

V^^%.^ {u, Vu) = V^[W„,pj^ (u, Vm)V^w^' + Wpi^ {u, VV;3«)] 

= Wuipi^ (w, Vu) V^^u-' + W^iukpi^ [u, Vu)V^u''V^u^ 
+ Wpip,^ {u, VV0u)Vl^iU^ + W„ipi^ [u, VVj3u)V^uK 

Then, integrating by parts and using Young's inequality, we have 
(3.22) 

VfiVaWpi {u, Vm)V/3 Au'dx= / VLWpi (m, Vm)V^^„m* 

° 7m2 ° 

^/ \W^ufdx-C [ |Vu|^(|Vu|^ + |V^u|2)da;. 

4 Jr2 iR2 

Note that |up = 1 implies 

u'V^ A + V;3w' Au' = -V^IVup. 



1R2 



> 
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By this identity, one can estimate the second term and the last term on the right- 
hand of (3.20) as follows: 

(3.23) 

/ V/3[Va {u''u%k {u, Vu)) - Vpk {u, Vu)w'^Vau']V/3 A dx 
= / ^Ib {u^Vp" {u, Vu)) u' Vff A u' dx 
+ [ V0u'Va{u''Vpk{u,Vu))Vi3 Au' dx 

<^ f \S/^u\^dx + C j IVupdV^up + IVu^da;. 

The other terms can be estimated easily in (3.20). Then it follows from (3.20)-(3.23) 

that 
(3.24) 

-I ^ r „ r 

73„,|2 



\^ f \Au\'dx + ^ I \VM'dx 

<C [ iVw^Vwl^ + \v\'^\V^u\^ + Q^^uf + |Vw|^)|Vw|2 dx. 
It follows from -u- Au= jVup, (3.18) and (3.24) that 
(3.25) ^ (y^^ + I V^l^) + I J^^ (IV^^P + iV^n dx 

<C [ ilvf + \Wuf) {IWvf + IV'up) dx. 
By the Gagliardo-Nirenberg-Sobolev inequality, one has 



C i\v\' + \\7u\')i\\7v\' + \V'u\')dx 

JR2 



<c( [ {\Wvf + \W\\'')dx] ' ( [ (IV^vp + IV^uHdx 



{{vl" + \Vu\'')dx 

/M2 



{\V\f + \V\f)dx+ (c [ (|t;|4 + |Vw|^)rfa; 

° JR2 V ^R* 

• (^J QVv\^ + \V^u\^)dx^ . 

This, together with (3.25), shows that for t e (0,T), 
(3.26) 

4/ (|V^;(a;,i)p + |VVa;,t)l')rfa;+^ / i\VM^,t)\^ + \^M^,t)f)dx 

"I JR2 O JR2 

<(c [ {\^u\* + \v\*)dx) I {\Vv{x,t)\^ + \yMx,t)\'')dx. 

\ ./M2 / Jr2 
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It follows from (3.9), (3.26), and Gronwall's inequality that for any s and t with 

T <s <t<T, 

(3.27) 

/ il'^vf + \V\\){x,t) dx 

< (eC/,7«2(|v«r+|.r)(x,t)d.d;^ (^j (|Vt;|2 + |v2w|2)(x,s)dx) 

< ("6^-1(1+™-^) / e{uo,vo)dx\ I (|V^;|2 + |V2w|2)(a;,s)rfa; 
Thanks to (3.8), (3.27), and the mean value theorem, we conclude that 

(3.28) sup / {\V'^u\ + \Wv\^){-,t)dx 

T<t<T Jm? 

<Ct-\1 + TR-^) E{uo,vo) eC-i(i+Tfl-^)E(«o,.o) 

for any t > 0. Then, by a similar proof as in Lemma 4, we can show that u belongs 

to C^/®(M^ X [t,T] for any r > 0. In the appendix below (Section 5), we can show 
that v) is regular for all t G (0, T]. □ 

Remark. Let u G V{0,T) and v G H{0,T) be a solution of (1.5)-(1.7) with initial 
values uq G H^{R^; S^), vq g -ff^(K^;M^) and divvo = 0. Assume that there are 
constants si and Rq > such that 

ess sup / |Vw(-,i)p + |w (-.t) p < El 

0<t<T,xeM2 Jbr{x) 

for any R G (0, Rq]. Then, for t G [0, T], we have 

(3.29) sup / {\V^u{x,t)\^ + \Vv{x,t)\^)dx 

0<t<T JR2 

< C3(l + Ti?-2) E{uo, vo) eCsi(i+™'')^^("o,-o)^ 

with Cs = C{\\uo\\h^ + \\vo\\m)- 

We are not able to prove the uniqueness of solutions to (1.5)-(1.9) for initial 
value in x as one in [Stl; Lemma 3.12]. However, we obtain 

Lemma 11. Let (mi, wi), (m2, ^2) £ V{i), T) x H{0,T) be two smooth solutions 
of (1.5)-(1.7) with smooth initial values {uo,vo) G Hl{m?;S'^) x jyi(R2;M2^ 
div Vq = 0. Then {ui,v{) = {u2,V2)- 

Proof. Following the proof of Proposition 15 in the Appendix, we can assume that 

|Vui| + IVU2I + + \V2\ < C 

for a constant C > 0. For simplicity, we set in (1.7) 

B{u, Vu) := - Wui {u, Vu) + W^k {u, Vu)u^u^ + W^k (u, Vu) VaW^^w' 
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It follows from (1.7) that 

(3.30) 



/ \{ui-U2)fdx+ [ {Wp^iul,Vul)-Wp^iu2,S7u2))S7^{u\-ui)dx 

= / {u\uWi (Ui, VUl) - ("2, VU2))V„(ui - 4) dx 

- / [-i3(wi, Viti) + B{u2-, Vu2) + {vi ■ Viti) - {v2 ■ V)w2] • {ui - U2) dx 

■■=h+l6- 

By Young's inequality, the last term on the right hand of the above identity can be 

estimated as 

h = - [-B{ui, Vui) + B{u2, VU2) + (f 1 • Vui) - {v2 ■ V)u2] • {ui - U2) dx 

<C \ui ~ U2\'^ + \vi - V2\'^ dx + ^ \V{ui - U2)\'^ dx. 

The difficult part is to estimate I^. Using an uniform open ball covering of M^, we 
can estimate only the local integral 

{u{u\V^j^{ui,yUi) - M^U2T^p3^(u2, Vm2))V(mi -U2)dx. 

Now we can think about in the equation (1.7) with = in a domain of M^. 
After a rotation TZ G 0(3), the integrand (1.2) has the following invariant property: 

w(7^u,7^Vu7^^) = w{u,Vu). 

Therefore, the system (1.5)-(1.7) is invariant for a rotation. Without loss of general- 
ity, we can assume that uo{xo) = (0, 0, 1). Since ui and U2 are uniformly continuous 
in {x, t) gR"^ X [0, t] for some r > 0, there exists a constant ro > such that for 
any {x,t) e Bro{xo) x [0,r] 



\ui{x,t) - uo{xo)\ < s, \u2ix,t) - ito(a;o)| < s. 



Then 



/ {u{u\Vpj^ (wi, Vui - Vu2)Va(ul - U2) dx 

J Brg{xo) 

<Ce j \V{ui-U2)\^dx + C j \W{u\-ul)\^dx. 



It follows from |w| = 1 that 
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Then an elementary calculation shows that 

|V(u? -ul)\< C\ui -U2\ + Ce\V{ui - ua)!- 
By a covering argument, we apply all above estimates to obtain 

h<Ce [ \V{ui-U2)\'^dx + C [ \{ui-U2)fdx. 

Therefore, choosing e sufficiently small yields 
(3.31) 

o37 / \{ui-U2)\'^dx + ^ I \V{ui-U2)fdx<C I \{ui - U2)f + \vi - V2\^ dx 

^ <M 7r2 I Jr2 J^2 

Using (1.5) and (1.6), one can obtain 

(3.32) 

\ii I \vi~V2\^dx+ [ \Vivi-V2)\'^dx 

^ Jr2 JR2 

<C \ {\vx - V2^ + |«i - U2^ + |V(ui - U2)r) dx + l- [ \V{vi - V2)\'^ dx. 

iR2 ^ Jr2 

Combining (3.31) with (3.32) gives 
(3.33) 

~f {C\u^-U2\^ + %i-V2f)dx<C [ {C\u,-U2f + %1-V2\^)dx. 
^ dt Jk2 4 J^2 4 

Integrating (3.33) in t and applying the Gronwall inequality, we conclude 

/ (C\u,-U2f + %i-V2\''){;t)dt<C [ {C\m-U2\'' + %1-V2\''){;0)dt = 0. 

JR2 4 Jg^2 4 

This proves our claim. □ 

4. Local existence and Proof of Theorem B 

In this section, we prove local existence of solutions of (1.5)-(1.7) and complete 
the proof of Theorem B. Recall the notation that V{T,t) denotes the space V{T,t) 
where is replaced by R^. 

Lemma 12. For a pair {uq, vq) e HKR^, S^) x L^iR^, R"^) with div vq = in in 
the sense of distribution, there is a local regular solution (ue, Vg) € V{0, T) x i?(0, T) 
of (1.10)-(1.12) with initial data (1.9) for some T > 0. 

Proof. Although Lin-Liu proved only the global existence of the solution to (1.10)- 
(1.12) with initial data (1.9) for the case of fci = ^2 = k^, their proofs still work for 
the local existence for the system (1.10)-(1.12). Thus we omit the details and refer 
readers to [LLl] and [LL2]. □ 
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Theorem 13. (Local existence) For a pair (uq. vq) e H^R'^.S'^) x L'^{M.^,R^) 
with div Wo = in K.^ in the sense of distribution, there is a local solution {u, v) € 
V{0,ti) X H{0,ti) of (1.5)-(1.7) with initial value {uq,vo) for some ti > 0. 

Proof. For any map mq S [M."^ , S"^) , one can approximate it by a sequence of 
smooth maps in H^(R^,5^). Without loss of generality, we assume that uq € 
H^(M.'^,S'^) and va e i?i(R2^]R3) ^j^-jj vq = in are smooth. Then thanks 
to Lemma 12, there is a local regular solution (w^, v^) € V x H oi (1.10)-(1.12) with 
initial data (1.9). 

For each pair {u,v), set 
ee{u,v) = W{u,Vu) + -^{l-\u\'^f + \v\'^, E{u,v)= I ee{u,v)dx. 

Then same calculations as for (3.1) give 
(4.1) 

E{Ue{-,t),Ve{-,t))+2 / / {\dtUe + {Ve-VUe)Ue\^ + \VVe\^fdxdt^E{uo,Vo)■ 
Jo 

By a similar analysis as in the proof of Lemma 8 and Lemma 7, one can show 

that there exist uniform positive constants Rq and £i, and a positive time = 
r(e, Rq. ci) such that the problem (1.10)-(1.12) with initial data (1.9) has a regular 
solution («£, «£) G V"(0,Te) x H{0,Ts) for each fixed £ > 0, and furthermore, it 
holds that 

(4.2) sup / \VUei;t)\'' + \Ve{;t)\'' + :^il-\u,{;t)\''fdx<ei 
0<t<T, JBr{xo) ^£ 

for any positive R< Rq. 

Next, we will show that there is a constant ti > 0, independently of e, such that 
> ti and the solutions (ue,We) is bounded in V{0,ti) x H{0,ti) uniformly in e. 

First, we claim that for all t & [0, min{l,T£}] 

(4.3) l<K{x,t)\<l. 

To verify (4.3), we re-scale the solution by 

u{x, t) = Ue{ex, e'^t), v{x, t) = evs{ex, e'^t), P{x, t) — e^Peiex, e^t). 
Then (u, v) solves the following approximate Ericksen-Leslie system 



(4.4) vi + {v ■ V)5' - Aff + V,,P = -V,, (V,,ti'= W^. {u, Vn)), 



(4.5) 



V-{; = 0, 
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(4.6) ui + {v- V){(* = V„ [Wpi^ {u, Vu)] - W^i {u, Vu) \uf) 
for z = 1, 2, 3, with initial data 

(4.7) v{x, 0) = vo{x), u{x, 0) = uo{x), Vx e IR^ 
where uo{x) = uo{ex) and vo{x) = evo{£x) satisfy 
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/ 



e{uo{x),vo{x)) dx = / e{uo{x),vo{x)) dx. 



M2 



The condition (4.2) becomes 

(4.8) esssup / \Vu{-,t){^ + U\^-\u{;t)\Y + \d\^dx<ei 

0<t<^,xeS?JBR{x) ^ 

£2 ' — 

for any R G (0,i?o]- While the basic energy identity (4.1) becomes 



(4.9) 



W{u,Vu) + -{l-\u\^f 



{■,t)dx 



Jm^ 



{\dt + \dtu + ii- V)un(-, t)dx dt = E{uo, vq) 



for all t e (0, |f). 

Without loss of generality, we assume < 1. Let r be the maximal time in 
[0, |f] such that 



(4.10) 



< \u{x,t)\ < 2. 



By (4.1) and similar arguments as for Lemma 8, one can derive from (4.4)-(4.6) 
that there exists a uniform constant Co such that 



(4.11) 

Co [ [ {\V^uf + \Vvf) dxdt 

Jo Jr2 

<-/ / Au-u{l-\u\'^)dxdt + C [ \Vuof + \vo\^ dx 
Jo Jm.2 

+ C [ [ \Vuf + \vfdxdt. 
Jo J-g? 

Integration by parts yields 

(4.12) -/ / Au ■u{l-\u\^)dxdt 

Jo Jm^ 

= 11 \'Slu\^{l~\u\^)dxdt-\ r f \V\u\^\'' dxdt 
Jo Jv? 2 Jq J-g2 
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By Young's inequality and using (4.1), (4.6) and (4.10), one can obtain 

/ / \Vu\'^{l-\uf)dxdt<Tj [ [ {l-\uffdxdt + c[ [ \Vu\^dxdt 
Jo Jo Jo Jr2 



JlR2 



{Vul'^dxdt + C 



I (IVt^ol^ 

JmP- 



+ \vo\^)dx 



10 Jv? 
for a small constant 

Combining (4.11)-(4.13) and choosing ei sufficiently small in (4.2) with Lemma 
2, we conclude that 

(iV^wp + IVCip) dxdt+ [ [ {l-\u\Y dxdt < C{1 + ^) I e{uo,vo)dx 



lo Jo Jm2 

for any R< Rq. 

It follows also from Lemma 2, (4.8), (4.9), and (4.14) that 

.2 ~ 



M2 



(4.15) 



Jm2 



{\Vu\^ + \dr)dxdt<Cei ( 1 + ;^ ) E{uo,vo) 



re 



for any R< Rq. 

Now following the calculation for (3.25), one can derive that for any t e (0,t), 
(4.16) 



<C I {\vf + \Vuf){\Vvf + \V\f){-,t)da 

JM2 



+ C 



M2 



V;3({(Xl-|wr))-V/j Au'dx 



+ C / |V{ir(-,t)da;. 

M2 



Note that 
(4.17) 



C 

= c 

< c 



I Wi3{u\l-\u\^))-Wi3 Au'dx 

JR2 

/ A{ff{l-\uf))- Affdx 

JR2 

/ \V^u{;t)\''dx + C [ \Vu{;t)\ 

JM2 JR2 



dx 



and 
(4.18) 



C [ \Vu{-,t)fdx = -C f \Wu{-,t)\'^u- Audx-C [ Va{\Vu\^)u ■ Vaudx 
Jm2 Jr2 J 

<?/ \^u{;t)fdx + C [ \Vu{;t)\''\V^u{;t)\dx. 
2 Jm2 Jm2 
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It follows from (4.16)-(4.18) that 

<C [ (|5|2 + |Vun(|Vi)|2 + |V2ti|2)da; + cf / {\V^u{;t)\^ + \^u{;t)\-')dx 

Using the Gagliardo-Nirenberg-Sobolev's inequality, one can get 
cf {\v\'' + \Wu\^){\Wi\'' + \V^u\^)dx 

<c( f + \\7u\^)dx)' ( [ {\Vvf + \V\\^)dx 



<c(^j ilVvf + \V^u\'')dxy (^J ilV^vf + \W^uf)dx^ 
Y {\v\^ + \Vu\^)dx^ 



<o/ il'^'W + \V'undx+{C i\vr + iVu^dx 



■ [ {\Vvf + \V^uf)dx. 

iR2 



This, together with (4.18), shows that 
(4.19) 



dt 

<(c 



^1 {\V'u\^ + \Vvr){.,t)dx+^ I {\V'u\' + \V'd\){;t)dx 



R2 



/ {\v\'' + \Vu\%,t)dx) [ {\V\f + \yvf){;t)dx + hit) 

with 

(4.20) h{t)=C f {\V^u\^ + \Vuf)i;t)dx. 

Jr^ 

It then follows from (4.14), (4.15), (4.19)-(4.20), and Gronwall's inequality that 

for all t e (0,t), 

(4.21) 

/ {\V^u\^ + \\/v\^){-,t)dx 

< ^CfX2(\v\' + \Vuf)(.l)<l^'il . f (|v2«o|2 + \Wvof)dx 

Jr2 

Jo 

< eC..(l+^)i^(«o,.o) (^ll^^lll^^ + WvoWl^ +J\{s)ds^ 
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Suppose that there is a xi e such that |{((a:i,t)| < 1/2 (or > |) 

with some t G [0,t]. It follows from (4.21) that u is C^-continuous uniformly in 
{x,t). Then, there is a constant C4 so that for x G Bi/4c^{xi) with < we 
have 



Then 

u 

.(^1) 



I [ il-\u{x,t)\ydx>l\B^m>e^ 

^ Jb 1 {X^) » 



which contradicts (4.8) for a sufficiently small ei. This shows that our claim (4.3) 
holds for all t G [0, 

Finally, we show that {ue,Ve) is bounded in F(0,min{T£, 1}) x H{0,mm{Ts, 1}) 
uniformly for any positive e < 4C4i?o- 

For any t < min(l,r£), it follows from (4.14) and (4.15) that 
(4.22) 

/ / {\V^Uef + \Wv'\''){x,t)dxdt+^f f {l-\Ue\''fdxdt 



(4.23) j J ^(IVuel^ + \v,\''){x,t)dxdt < C£i (^1 + ) E{uo,vo). 



Let (f be the cut-off function as in the proof of Lemma 9. Then by a similar 
analysis as in (3.14)-(3.17) and using (4.22)-(4.23), one can get 

(4.24) 

/ e{ue{-,t),v,{-,t))ip'^dx+ I I {\Vve\'^ + hdtU,+v,-Vue\'^)ip'^dxdt 

Jr^ Jo Jr'^ ^ 

< j ^e{uQ,VQ)ip^ dx + C^—(^ + E{uo,vo) 

+ fl ^{^-W?f\Ve-^{f)\dxdt. 

Jo Jr2 ze 
On the other hand, 

t r 1 



(4.25) / / il-\u,\'y\v,-\7iip')\dxdt 

Jo iR2 ^£ 
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where one has used (4.1) and (4.22). 

Hence, 
(4.26) 

/ ee{us{x,t),Vs{x,t))dx < / e{uo,vo)dx + C^— [1 + ^] E{uo,vo) 

for any R < Rq. First, choosing i?i > so that 
(4.27) / e{uo,vo)dx <^ 

for all xo e M^. Then, set 

£1 i?i 



(4.28) ti=m:m{Ri, 



4C(£i +£;(uo,t^o)) 



Then for < < ii, 

/ ee{Ue{x, t), Ve{x, t))dx < £i 

JBr{xo) 

for all X(j e and R < Ri. Consequently, we have shown that there is a uni- 
form ti < min{Tg, 1} such that {ue,Ve) is bounded in V{0,ti) x H{0,ti) with ti 
independent of s. Letting £ ^ 0, we can prove the local existence of solution 
{u,v) G V{0,ti) X H{0,ti) with initial data (1.9). □ 

Now we complete the proof of Theorem B. 

Proof of Theorem B. By Theorem 13, there is a local solution {u,v) G V{0,ti) x 
H{0,ti) of (1.5)-(1.7) in x [Q^ti] with initial conditions (1.9) for some h > 0. 
By Lemmas 10-11, the solution can be extended in [0,Ti) for a maximal times Ti 
such that at Ti, there is at least a singular point xj G such that 

limsup / e{u,v){-,t) dx > sq 
tyTi Jbr{x]) 

for any R < Rq for some i?o > and So > 0. It is easy to see the solution 
[u.v) £ V X H ]s regular for all t G {Q,Ti). Then there exists a sequence of 
{t„} such that the sequence (u(t„), f;(t„)) converges weakly to {u{T\),v{Tx)) in 
H\M?-S^) X L2(R2.]g3^ satisfying 

/ e{u{Ti)^v{Ti)) dx < I e{uo,vo) dx — Sq, div i;(ri) = 0. 

Using the energy identity, there is a finite number of singular times {Ti}^-^ in 
Theorem B. □ 
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5. Appendix: The liquid crystal flow and regularity issue 



In this section, we formulate the liquid crystal heat flow and discuss C^'"- 
regularity issues for solutions of the liquid crystal flow (1.4) and the system (1.5)- 
(1.7). 

The liquid crystal equilibrium system in a form of vectors and tensors was derived 
by Hardt, Lin and Kinderlehler in [HLK] using the Lagrange multiplier method, 
but we need a precise form of (1.3) in coordinates. 

Let (I) he a smooth functional in Co°(0,R^). We consider a variation 



utix) 



U + t(j> 



U + t(f) 



\U{X) + t<l){x)\ (l + 2to-^ + t2^2)l/2 



and compute 



dut 
~dt 



{u + t<p){u-<p + t\(j)\^ 



(l + 2to-.^ + t2</,2)l/2- 

To derive the Euler-Lagrange equations, we compute 



^ jjV{uuVut)dx 



0. 



t=o 



This implies 



dx = 0, 



t=o 



where Wp^{u,p) = ^ and W^. = ^. Note 



du\ 
dt 



= 4>^ - u^{u ■ 



dVaUl 



t=0 

We conclude that 



dt 



t=0 



(5.1) [w^,{u,Vu)[cl^-u^{u-cj>)] 
Jn 



+ Wp^Ju,Vu) [Vc.(/)* - \/cu\u- (j)) -u'\7a{u■ 



dx=0 



for any (j) e C(f'(il, R'^). Therefore, we call that u € H^{Q,S'^) is a weak solution 
to the liquid crystal system if u satisfles 

- Va [Wp.^ (u, Vu) - u'^u'Wpk^ {u, Vli)] + Wu^ {u, Vu) 

- W^k (u, \7u)u''u' - Wpk^ {u, Vu) Vaw'^M* - Wpk {u, yu)u^Vo,vj- = 

in the sense of distribution. Note \u\^ = 1, then u'Vu* = 0. This system is the 
exact form of (1.3). 

Then, the liquid crystal flow can be formulated as in (1.4), i.e.. 



'dt 



=Va [Wpi^ {u, Vu) - u'^uWpk {u, Vu)] - W„i (u, Vu) 
+ W„fe {u, Vu)u^u'' + W k {u, Vu)Vc<u''u' + V k {u, \/u)u'"Vau\ 
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Next, we will prove that a Holder continuous solution of (1.4) belongs to C^'" 
for some a with < a < 1. For any point zq = {xq, to) S O x [0, M) and any number 
i? > 0, we use standard notations: 

Bixa,R) = {x e : \x - xo\ < R}, Q{zq,R) ^ B{xq,R) x {to-R^M), 
Sr{zo) = B{xo, R) X {to - i?'} n dB{xo, R) x {to - R^, to). 

Proposition 14. Let fl be a domain in M"^ with smooth boundary dfl. Let u be a 
weak solution of (1-4) and Holder continuous in Q, x [0,T). Then, Vu is (locally) 
Holder continuous with the same exponent in Q x [0,T). 

Proof. Assume that u(x,t) is Holder continuous with exponent /3, < /? < 1. 
Let {xa,ta) e Q x {0,T) with QiB.,{zQ) C x (O.T) for some Ra > 0. Note 
u{xo, to) = e G S*^. After a rotation, we can assume that e = (0, 0, 1). 

It follows from |u| = 1 and Cauchy's inequality that 

(5.2) |u3|2 |Vu3|2 < _ |u3|2)|yy|2 < 2\u - u{xo,to)\\Vu\'' . 

Denote 

P=(pL)3x2. 
Using the structure of W{u,p), we can write 

Wp{u, Vu) = Wp{u, Vu\ Vu^) + f{u, Vu^), 

where \f{u,Vu^)\ < C|Vu^|. 

Let V = {v^,v^) be the solution of the Cauchy-Dirichlet problem 

(5.3) vl = V„ [Wp.^ (e, Wv\Vv^)] in Qn{zo) 
v' = on Sr{zo). 

for i = 1,2. Since (5.3) is a parabolic system with constant coefScicnts, it follows 
from Proposition 1.2 in [GS; Proposition 1.2] that for all p < -R < -Ro 



and 




Set w = u — V. Then for all p < R,we have 

(5.4) / \S/ufdz<c(^Y [ \Vu\^dz + C [ \Vw\^dz + C [ \Vu^\^ dz 
JQp JQr JQr JQr 



34 



M.-C. HONG AND Z. XIN 



and 
(5.5) 

/ \'^u-{Vu)p\'^dz<c(^y [ \Vu-{Vu)Rfdz + C [ I'^wfdz 
JQ, JQr JQr 

+ C [ \Vu^\'^dz. 
JQr 

Note that u is /3-H61der continuous in f2 x [0,T) and u{xo,tQ) = (0,0, 1). 

Although there is no maximum principle for the parabolic system (5.3) with 
constant coefficients, Giaquinta-Struwe in [GS; page 445] obtained that 

sup \v — u{xo,tQ)\ < Csup \u — u{xq, to) I 
Qr Qr 

with a constant C independent of R and u. This implies 

(5.6) \w\ < \u-u{xo,to)\ + \v-u{xo,to)\ < CR^ . 

Multiplying the difference between (5.3) and (1.4) by {i = 1,2) and integrating 
over Qr lead to 

(5.7) 

2 

/ \w\^{-,tQ)dx+ / ^S/ocW^Wp^^{e,S/w)dx 

Jbr JQr 

<[ y^\Vaw'\\Wp^{e,yv)-Wp^{u,Vu)\dx + C I \Vu^\\Vw\dx 
JQr i=i " " JQr 

+ f ^^Va^i;*wV=ypj(w,Vw) + C /" \w\\Vu\^dx. 
■'Qr i=i fe=i " ''Qr 

Since u is /3-H61der continuous and u(xo,io) = (0,0,1), we have < CR^ for 
i = 1,2. Applying Young's inequality and (5.2) yields 

(5.8) / \Ww\^dz<CR'^ [ \Vu\^dz. 
JQr JQr 

It follows that for all p < R, 

(5.9) / \Wu\''dz<c(^Y [ \Vu\^dz + CR^ j \Vu\^ dz. 
jQp JQr JQr 

We claim the following Cacciopoli's inequality 

(5.10) / |Vuprf2<C^ / \u-U2R\^dz<CR^+^^. 

JQ{zo,R) JQ{zo,2R) 

for any zq €0, x (0, oo) and R < Rq, where U2r is the average of u in Q2R{xo,to). 
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Next, we prove this claim. Let ^ be a cut-ofF function in C^{B2r{xq)) with 
0<C<l,C = lin Br{xo) and |VC| < Let r e C°°(R,R) be a function depends 
only on t with < r < 1, r = 1 on [io — -R^, <o] and r = on (— oo, to — 4:E?) and 
\dtT\ < C/B?. 

Testing (1.4) with (p = {u^ - '"2_R)C^'r^-f(-oo,to) for i = 1,2, where I{-ao,to) is the 
characteristic function of (— oo, ^0), we have 

2 

/ \U{; to) - U2R\^eT\to) dx+ [ V W^, , ^aU'VaU'^^T^ dz 

<2 f [Wp^ {u, Vu) - u''u%k {u, Vu)] Va^(u' - dz 

JQ2r(zo) 

+ C j \Vuf\u-U2R\fT^dz + 2 f \u-U2R\''^^TdtTdz 

+ f u''u%kJu,Vu)Va{u' -uiR)S,^T^dz + C f IVw^l^fr^dz. 

JQ2r{zo) " JQ2r(zo) 

Since u is /3-H61der continuous and u{xq, to) = (0, 0, 1), u{x, t) — U2r can be chosen 
sufficiently small when i?o is small and + \v?\ is also small. We need to deal 
with the above last term. By (5.2), the term |Vu^p is also good. By Young's 
inequality, the claim (5.10) is proved. 

Using (5.9) and (5.10), a standard iteration (cf. [G], Chapter III, Lemma 2.1) 
yields that for all p < Ro, one has 

(5.11) / IVwl'dz < Cp3+3/3, 

where C depends on Rq. An iteration by (5.9) and (5.10) yields that for any a <1, 

[ \Vufdz<Cp^+^^. 
JQp 

Using (5.2) and (5.8) yields 

/ \Vu-{Vu)pfdz<c(^'' [ \Vu-{Vu)R\^dz + CR^ j \Vu\^ dz 

JQp JQr JQr 

<c(^y f \Vu - {Vu)Rf dz + Ci?3+2-+/3. 

jQn 



Choose a sufficiently close to 1 so that 2a + ^ > 2. Then, for all p < -j, we have 

/ \Wu-{Wu)pfdz<Cp^+^''^ 
JQp 

for some ai with < cti < 1. This implies Vu € Cf^^^ and then Vu € C^'^ (cf 
[GS]). □ 
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Proposition 15. Let {u,v) be a weak solution of (1.5)-(1.7) in x [0,r] and 
assume that u is Holder continuous in x [0,T). Let r be any positive constant. 
Fort G [t,T], we have 

f \W^u{x,t)\'' + \Wv{x,t)\^dx < Ct-^{1 + TR-^). 

Then, {u,v) is smooth in x (0,r). 

Proof. By Sobolev's embedding Theorem, we have 

\Vu{x,t)\P + \v{x,t)\Pdx 



'Si(xo) 

for any p > 1 and for xq G and t > t. By a similar way to one in Lemma 5, we 
can show that u is Holder continuous in x [r, T] . 

To get the higher order regularity, we rewrite (1.7) as 

(5.12) ui - Va [Wpi^ {U, Vu)] = -u'^u'Va [Vpk {u, Vu)] - {v ■ V)u' + B{u, Vu), 
where B{u, Vu) is given by 

B{u, Vu) = - Wu^ (u, V?i) + {u, Vu)u^u'^ + Wpk {u, S/u)Vau''u'' 

+ Vpfc (m, \/u)u''VaU^ - Va [u'^U^] Vpk {u, Vu). 

Since W{u,p) is quadratic and convex in p, we can write 

WpiJu,Vu) = ai^^{u)Vau\ 

Since u is uniformly Holder continuous, the left-hand term of (5.12) is a parabolic 

operator. Let £_{x) be a cut-off function in Bji(xo) and let r G C°°(R,R) be a 
function depends only on t with < r < 1, t = 1 on [to — j-R^, io] and r = on 
(-00, to - R^) and \dtT\ < C/R^. Set (j) = r^. Multiplying (5.12) by <f>, we have 

(5.13) {u<l>)l - [a%{v)V^{u^^)\ - u'ct>t 

= -u'^u'Va [Vpk {u, Vu)] cl>-[{v- V)u* + B{u, Vu)]0. 

By the assumption, we have 

{vV)uGLPiQnixo)), \Vu\'' e LPiQnixo)) Vp > 1 

But the first term on the righthand of (5.13) is not a 'good' term, which need more 
analysis. Using the fact that |u| = 1, we have 

u^VlpU^ = -iypu ■ Vu + u^VipU^ + u^V'ipU^), u^uf = -{u'^ul + u^ul). 

Without loss of generality, we regard the solution in R^. By a rotation, we assume 

u{xqM) = (0,0,1). 
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Since u is Holder continuous, there exists a small R such that 

\u{x, t) — u{xo, to)I < e 
for a sufficiently small constant e > 0. Therefore 

IV'w^i < c-IVul^ + 2£(|VV| + IV^u^l) 

Apply the classical -L^-estimate of parabolic systems (c.f. [Ei], [LSU]) to (5.13) 
for i = 1,2, we have 

WutnLHQnicco)) + I|V'M)||lp(Q«(.o)) < CUV\^L. + Cs\\cl>W\\\ l 

+ <^(l|w||L2f(Q„(a;o)) + \\v\\l^p{Qi,{xo)) + 1)> 

where u = (v},v?'). Choosing e sufficiently small, we obtain 

\\ut(t>\\Lv{Q^{xo)) + \\^'^{'^4>)\\lv{Qr{xo)) < C. 

To estimate v in (1.5), it follows from Holder's inequality that 

[ \{v V)v\P dx < i [ \Vv\^dxdt] if \v\^dxdt] 

for any p with 3 < p < 4. By the L^-estimate of Stoke's operator (e.g. [So]), Vt and 
V^u are in for 3 < p < 4. This implies that v is Holder continuous. 

Differentiating in xi in (5.12), we have 

= -u'^u'Woc \ypi {u, VVx, w)] + u#V^M + Vt;#Vu + Vu#V^u. 

By applying the i^-theory, a similar argument yields that Vw is uniformly contin- 
uous. Then, a standard bootstrap method implies that {u, v) are smooth. □ 
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